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Typical examples of Nih) data from a series of NBS publications on ionospheric electron 
densities are described briefly as an introduction to data available from a one-year's program 
beginning during the International Geophysical Cooperation, 1959. The entire body of 
data is then illustrated in compact form, and discussed phenomenologically. In a synthesis 
and interpretation, it is concluded that diurnal, seasonal, and latitudinal temperature 
variations in the F region may explain many features of the quiet-day behavior of that 
region. A corpuscular component of heating at mid and high latitudes is suggested as 
accounting for the seasonal anomaly in the daytime F region; the anomaly is assigned to 
the summer season rather than to winter, on the basis of evidence given. The nighttime 
electron density variations are found to be explainable by the loss rate at the equilibrium 
height to which the layer drifts under the influence of diffusion. An appendix discusses 
the day-to-day variability of the data. 



1. Introduction 



By the end of tlie IGY, the number of ionospheric 
soundings reduced to electron density profiles was 
comparatively snudl, measm*ed in terms of the very 
large vertical soundings program then underway. 
The original difficulties stemmed from the consider- 
able labor necessary for manual computation of 
single profiles, but by the time of the IGY, the 
availability of electronic computers had largely 
eliminated this. The remaining limitation was the 
rather formidable amount of scaled data necessary 
for each computation. The IGY electron density 
profile surveys conducted by Schmerling [1957] and 
Thomas and Vickers [Thomas, 1959] for selected 
stations on Intei-national Quiet or Regular World 
Days were limited in their scope by the large amount 
of data preparation necessary. 

When the National Bureau of Standards embarked 
upon a similar program early in the IGC (1959), 
part of the effort was devoted to developing tech- 
niques by which the valuable manpower resources 
of the network of U.S. and associated vertical 
soundings stations — ^already the mainstay of the 
conventional ionospheric soundings data published 
by NBS — could be applied to numerical reduction 
of their ionograms for N{h) analyses. Simple 
methods were developed for this purpose, permitting 
the station scientist to derive and tabulate data for 
subsequent N(h) analysis by the central laboratory's 
computer. These methods, applied to the hourly 
ionograms of a day's observations, require about 
the same effort as the preparation of a quarter- 
hourly /-plot. 

Since May 1959, the NBS has conducted an 
hourly electron density profile survey, using data 
provided by the ionospheric stations, from a current 



total of 11 stations. These stations, their affiliation, 
and their initial date of participation in the program 
are given in table 1. 

Our first objective in the study of these data has 
been to classify and describe the mean quiet geo- 
graphical, temporal, and Iveight structure of the 
northern mid-latitude ionosphere. For this purpose, 
five of these stations (Newfoundland, Ft. Monmouth, 
White Sands, Grand Bahama Island, and Puerto 
Rico) are close enough together geographically to 
permit the delineation of the structure of the mean 
quiet ionosphere between geographic latitudes of 
15° to 50 °N, and geomagnetic latitudes of 30 °N 
to 59 °N, with some confidence. The results from 
these stations for one year of data (March 1959- 
April 1960) are being issued by tlie NBS in a series 
of Technical Notes [Wright et al., 1959-1961]. It 
is the purpose of this paper to review this considerable 
volume of data, to portray it in a different summary 
form, and to discuss certain other results which 
assist in the interpretation of the whole. Finally 
there will be given a synthesis and interpret {it ion 
of the data in the light of contemporary theory. 







Table 1 








station 


Affiliation 


Latitude 


Longitude 


Mag, dip 


Began 

mil) 


Puerto Rico 

Grand Bahama Is.-_ 

Fort Monmouth 

White Sands, N.M., 

St. Johns, Nfd 

Adak, Alaska 

Okinawa, Ryukus.- 
Thule, Greenland. ._ 

Huancayo,Peru 

Talara, Peru 

Baguio, Philippmes. 


NBS 

USASigC. 
USA SigC _ 
USASigC. 
DRTE__._ 
USASigC. 
USASigC. 
USASigC. 

NBS 

NBS 

NBS 


18°30' N 
26°40' N 
40°15' N 
32°24' N 
47°33' N 
51°54' N 
26°30' N 
76°3r N 
12°03' S 
4°34' S 
16°25' N 


67°12' W 
78°22' W 
74°()1' W 

10()°52' W 
52°4()' W 

176°39' W 
128° W 
f)8°50' W 
75°20' W 
81°15' W 
120°36' E 


5L5° N 

59.5° N 

5L7° N 

60° N 

72° N 

63° N 

36.5° N 

86.2° N 

0.5° N 

13° N 

17° N 


Jan. 1959 
Feb. 1959 
Feb. 1959 
Mar. 1959 
June 1959 
June 1959 
June 1959 
July 1959 
Jan. 1960 
Jan. 1960 
Feb. 1960 
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2. Calculation of the Electron Density Data 

For the data discussed here, the well known 
matrix method of Bndden [1955] has been used to 
derive the electron density profiles. It is unneces- 
sary to review here the details of the primary virtual 
height-to-true height conversion, since in most 
respects the procedure has corresponded exactly to 
Budden's. However, it should be observed that 
this method has several shortcomings which lead to 
errors in the N{h) data. 

The shortcomings are nearly as well known as the 
method itself and do not need detailed discussion 
here; it suffices to itemize them: (a) The electron 
density is assumed to be a monotonic function of 
height ; if a valley exists there will be a height error 
which diminishes from a value equal to the valley 
width just above the valley, to smaller values at 
greater frequencies, (b) The shape of the profile is 
somewhat inaccurately determined in regions of 
strong curvature, because of a convenient but 
inappropriate assumption made about the variation 
of electron density within small intervals [Paul, 
I960], (c) Because the ionogram observations them- 
selves do not ordinarily contain information about 
electron densities below about 10^/cm^, there is 
usually an error at night due to neglect of retarda- 
tion in this ionization. 

While one or more of these may be serious for 
individual profiles, and while their resultant may give 
mean profile parameters which contain net errors of 
perhaps 10 percent, these errors are small compared 
with the total variations discussed here and should 
not invalidate our conclusions. 

The result of the primary matrix multiplication is 
a table of true heights at particular plasma frequen- 
cies. It is more convenient to derive from these 
data certain other quantities amenable to convenient 
ph3^sical interpretation. Our primary true height 
data has therefore been put through a secondary 
calculation process in which the following parameters 
are derived : 

2.1. Electron Density at Fixed Heights 

The plasma frequencies (Jn-) are converted to 
electron densities by the relation A^(electrons/cm^) 
= 12,400/^(Mc/s)^. A linear interpolation among the 
corresponding true heights then provides the electron 
density at 10 km height intervals throughout the 
observable portion of the profile. 

2.2. Height and Characteristic Thickness at the F2 

Peak 

lonograms contain no direct information from a 
layer peak itself, since the virtual height of a radio- 
frequency penetrating just to this level is immeasur- 
able. Nevertheless, this level is of unique interest, 
and it is essential to devise means for its description. 
A practical method is to fit the portion of the true 
height profile near the peak with a suitable curve, 
and to determine the parameters of the peak from 



this curve. The parabola is the simplest curve for 
this purpose, and has the additional merit of closely 
approximating the peak of a ^ ^Chapman'' distribu- 
tion. The parabola is given by 
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where /imax is the height of the layer peak, iVmax 
the peak electron density corresponding to the critical 
frequency, and Ym is a parameter characteristic of 
the thickness of the layer. This curve has been fit 
to the highest portion of our true height curves, using 
the measured critical frequency and two true heights: 
the highest and the fourth from highest. The ac- 
curacy with which the parabola may be fit to these 
data depends slightly upon the spacing of the two 
heights and rather critically upon the accuracy of de- 
termination of A^max. In the data described here, 
various checks have been applied to eliminate ex- 
treme errors due to the latter cause. The frequency 
spacing of the two true heights used for fitting the 
parabola is a compromise between a narrow spacing 
(sensitive to small relative errors in the data) and a 
wide spacing (sensitive to real departures from the 
parabola). In practice, these points are separated 
by about 0.8 Mc/s. 

In addition to the height of the peak, the quantity 
Ym/2, which we call "Scaf\^ is determined. This 
is the quarter-thickness of the parabola, and is taken 
as a measure of the scale height of the atomic species 
at the level Am ax i^2, in what follows. The reasons 
for this interpretation are given briefly in the follow- 
ing section. 

2.3 Extrapolation of Profiles Above hmax F2 

For sunspot maximum conditions it was previously 
shown [Wright, 1960] that an electron distribution 
above the F2 peak of the form 



1 The terms Scat, Shmax, -S/iinf have evolved in our N(h) system as short pro- 
nounceable names for Scale Height, "/"to hmax" and "/to hoo", respectively. 



agreed fairly well with available rocket and other [ 
data, for a neutral particle scale height of H=100 
km. This is, of course, the well-known equation 
derived by Chapman for the equilibrium distribution 
due to electron production and recombination-type 
loss, but it is also the equilibrium form theoretically 
expected of a layer under the combined influences 
of diffusion and attachment-like loss, as has been 
shown by Hirono [1955], among others. Since these 
latter processes are thought to be effective in deter- 
mining the F2 peak [Ratcliffe et al., 1956], the use 
of eq (2) for extrapolation above the F2 peak has 
some theoretical justification. Since eq (2) approxi- 
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mates to the parabola of eq (1) with Ym = 2H, the 
interpretation of Scat in terms of the scale height 
receives a similar (l(*o;rcc of theoretical justification. 
Each of the hourly profiles used in this study has 
been extrapolated aboYe hiniixF2 by eq (2), with 
H=100 km. The total electron content, which we 
term Shm(, is given by adding to the subpeak con- 
tent Shmax (vide infra), the quantity 2.82^ A^max, 
as may be shown easily by integration of (2) between 
limits Amax and infinit3^ 

2.4. Subpeak Electron Content 

The total number of electrons in a unit column 
extending from the lowest observable ionization (usu- 
ally about 10*/cm^) to /?,max, is termed Shmax. It 
has been obtained in the course of our secondary 
calculations by numerical integration of the profiles 
between these two limits. 

2.5. Mean Electron Density Profile Data 

The systematic calculation of hourl}^ profiles pro- 
vides a great quantity of data containing, among 
other inhomogeneities, various degrees of ionospheric 
disturbance. Our first objective has been to derive 
mean quiet conditions from these data. This has 
been done by eliminating profiles at those hours 
for which the magnetic character figure Kp exceeds 
4; from the remaining data, the montldy mean val- 
ues of electron density at 10 km height intervals, 
and similar averages of the other special quantities 
are obtained for each hour. Generally, data from 
about 20 profiles comprise such a mean. 

At the same time, the standard deviation of the 
data entering each mean is obtaincMl. For some of 
the parameters, a more useful measure of the vari- 
ability is the relative standard deviation or per- 
centage variability of the quantity. This is obtained 
by dividing the standard deviation by the mean. 
From the program described here, the behavior of 
ahmax, (aScat)/Scat, (aShnnix) /Shumx, and (o-iVmax)/ 
A^max is discussed in an appendix to this paper. 



Typical Examples of TN 40- 
Representations 



-Series 



Tlie ]nean values of the quiet-day electron den- 
sity profiles and their derived parameters, for each 
month of the year March 1959-February 1960, are 
portraved in a variety of representations in the 
NBS Technical Note series 40-1, 2, 3, ... . etc. 
Typical diagrams from this series will first be illus- 
trated and discussed here, as an introduction to the 
parameters and as a brief review of the variety of 
information available. While conmients are offered 
regarding the significance and interpretation of spe- 
cific phenomena evident in these diagrams, no at- 
tempt is made here to discuss the generality of the 
phenomena; rather, the aim is to point out several 
of the special properties of each form of representa- 
tion. Throughout, it should be borne in mind that 
our data and conclusions pertain to a period near 
the maximum of the solar cycle. In the following 



sections the entire body of data will then be repre- 
sented in a series of month-by-hour diagrams (for 
a year running from May 1959 through April 1960), 
which summarize the diurnal and seasonal variations 
in the northern mid-latitude ionosphere. 

3.1. Vertical Cross Sections 

Figure 1 (from NBS Tech. Note 40-6) illustrates 
a vertical cross section of the ionosphere, nominally 
above the 75° W meridian between geographic lati- 
tudes of 15° and 50° N, for 1700 75° W time, Au- 
gust 1959. Contours of plasma frequency /at in Mc/s, 
related to electron density N by 12,400/^— A^ elec- 
tron/cm^, represent the true heights of reflection of 
vertically incident radio waves. The height of max- 
imum electron density is represented by the dashed 
line, and electron densities above this level are the 
result of an extrapolation according to the model 
discussed in section 2.3. Lines of the geomagnetic 
field would be very nearly vertical within this dia- 
gram, and are therefore omitted. Note tliat the 
vertical scale is expanded relative to the horizontal 
by a factor of about 5.5. 

Near the level of maxinmni density a latitudinal 
gradient of about 0.0015 Mc/s per kilometer or 0.16 
Mc/s per degree of latitude may be seen, and this is 
a fairly typical maximum value at this latitude. It 
becomes rapidly smaller at heights below the jnaxi- 
nmm, and not so rapidly smaller at heights above the 
maximum. It is interesting to note that the latitu- 



AUGUST 1959 
1700 LMT 




45" 40' 35° W 25" 20" 

NORTH GEOGRAPHIC LATITUDE 



Figure 1. Vertical cross section along 75° W meridian. 
Contours of plasma frequency in Mc/s. 

Dasiicd line rcprosonts height of F2 peak. Contours above this level obtained 
from Chapman model with H=100 km. 
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dinal gradient considerably exceeds the longitudinal 
sunrise gradient. A typical value of the sunrise 
gradient is 0.0007 Mc/s per kilometer or 0.07 Mc/s 
per degree of longitude near the F2 layer maximum. 

3.2. Electron Density Maps in Latitude Versus Time 

Contours of electron density at constant heights 
of 200 and 300 km are shown in figures 2 and 3, for 
the months June 1959 and January 1960, respectively 
In these diagrams, the contours are given at intervals 
of 10^ electrons/cm^. In addition to illustrating the 
diurnal variation of electron density at a fixed 
height versus latitude, these diagrams may also be 
considered as representing the mean longitudinal 
variations when it is noon over the 75 °W meridian. 
For this purpose, the latitude scale is expanded by 
a factor of 10.76 relative to the ^^ongitude^' (time) 
scale. 

At 200 km (fig. 2), the almost exclusively solar 
control of the ionization densities is obvious: the 
electron densities rise rapidly at sunrise from values 
less than lOVcm^ (the minimum detectable with 
ordinary observing devices) and decrease below this 
limit a little more slowly at sunset. A small latitu- 
dinal gradient may be seen corresponding to the 
greater solar zenith angles nearer the equator. Note, 
however, that at this season one of the anomalies of 
the ionosphere at higher altitudes is reversed at 200 
km: the maximum electron density is reached hejore 
noon, at least at the higher latitudes. This diurnal 
asymmetry has been noted by Croom et al. [1959]. 
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Figure 2. 
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ELECTRON DENSITY AT 300 KILOMETERS 
JANUARY I960 




Latitude versus local-time isoionic map at 200 km. 
Contours in electrons per cm3X10-\ 



CONTOURS = el/cm^xlO"^ 



Figure 3. Latitude versus local-time isoionic map at 300 km. 

Contours in electrons per cmsxiQ-^. 

At 300 km (fig. 3) strong solar control is also 
evident, but it is clear that other factors are rela- 
tively more important than at lower altitudes. The 
initial increase of electron density at sunrise is nearly 
the same at all latitudes, although the later midmorn- 
ing values are considerably greater at low latitudes. 
There is also a peak above 45° N, which is partially 
explained by the lower altitude of the entire layer 
there. The curious irregularity in the midday con- 
tours, which appears to ^^inove'^ to later hours at 
lower latitudes is explained by the latitude/time 
variations of the height of maximum density: 
Amaxi^2 is nearh^ 300 km and varies with time above 
and below 300 km in the vicinity of this irregularity. 

3.3. Latitude Versus Time Variation of i2maxF2 

The height of the F2 peak is illustrated for March 
1959, in figure 4. This is one of the simplest of 
similar diagrams for other months. Its most obvi- 
ous feature is that the F2 peak falls at sunrise from 
about 380 km to about 300 km uniformily at all 
latitudes. The daytime variation is characterized 
by a gradual rise at all latitudes, but this proceeds 
more rapidly, and reaches a higher daytime altitude, 
at the lower latitudes. It continues to increase until 
sunrise at middle and higher latitudes, but decreases 
somewhat before midnight at low latitudes. 

3.4. Latitude Versus Time Variation of F2 Layer 
Thickness 

The ^^thickness^' of the F region is characterized 
by our quantity ^^Scat/^ the quarter- thickness of a 



300 



N^SO 



HEIGHT OF MAXIMUM ELECTRON DENSITY 

MARCH 1959 




—I — I I i_ 
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CONTOURS =km 



Figure 4. Latitude versus local-time map of the height of the 

F2 peak. 

Contours in kiloniotcrs. 



p^vrabola fit to the F2 peak. Its lalitiuiiiial and time 
variations arc illustrated in ligure 5 for the month 
of July 1959. This parameter is of special interest 
since it is probably an approximate measure of the 
neutral particle scale height — and lience the tempera- 
ture — at the F2 peak. It is at once clear that the 
thickness is strongly under solar control, increasing 
at sunrise from a nighttime value of about 50 km to 
about 70 km during da^^time. In a significant 
departure from solar control, however, the higher 
latitudes have a thicker daytime F region than the 
lower latitudes. 

3.5. Latitude Versus Time Variations of Sub-Peak 
Electron Content 

The number of electrons in a unit column below the 
F2 peak is designated ASAmax; its latitudinal versus 
time variations for May 1959 are illustrated in figure 
6. The behavior of the electron content is extremely 
simple, closely following the solar zenith angle. 
Nevertheless, certain anomalies are evident. At 
low latitudes, the maximum value of S/^max follows 
noon by l}^-2 hours, rather in accord with simple 
theory, while in middle latitudes, the maximum is 
distinctly at noon. It appears that at higher lati- 
tudes the maximum again shifts to post-noon. Also 
the latitudinal midday gradient in /SAmax is very 
much larger than would be expected from the lati- 
tudinal change in midday solar zenith angle. 



QUARTER THICKNESS OF F-REGION PEAK 
SCAT 

JULY 1959 




CONTOURS = km 



Figure 5. Latitude versus local-time map of the quarter- 
thickness of the F2 peak. 
Contours in kilometers. 



ELECTRON DENSITY INTEGRATED TO 
HEIGHT OF MAXIMUM ELECTRON DENSITY 

MAY 1359 




12 
LMT 



18 00 

CONTOURS = el/cm^ col xlO~'^ 



Figure 6. Latitude versus local-time map of the subpeak 

electron content. 

Contours in electrons per cm2 columnXlO-12. 
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3.6. Diurnal Variation of Electron Density at Fixed 
Heights 

The diurnal variation of electron density at a given 
heiglit is at first sight the most natural representa- 
tion of ionospheric variations. An example is 
shown in figure 7, for Puerto Rico, April 1959. 
When the height of maximum falls below the fixed 
heights represented here, the corresponding curves 
are dashed. Such electron densities are the result 
of tlie topside model discussed above. The envelope 
of these curves corresponds to the variation of the 
maximum density. 

It is interesting to note the progressive departure 
from simple solar control at successively greater 
altitudes. While the sunrise period (and until 
1000) is characterized by a rapid increase of electron 
density at all heights, there then occurs a general 
decrease, most marked at heights above 250 km. 
That at least part of this behavior might be attribut- 
able to movement of the layer, is suggested by the 
post-sunset increases in electron density which occur 
at successively later times at lower altitudes. 



N(t) PUERTO RICO APRIL 1959 
FIXED HEIGHTS IN KILOMETERS 




Figure 7. Diurnal variation of electron density at fixed heights. 
Contours in kilometers. 



4. Diurnal and Seasonal Variations 

A great deal of the basic information describing 
the F region of the ionosphere is contained in the 
parameters hmsixF, Scat, A^max, and Shmo^x, which 
measure the height, thickness, and electron density 
at the peak, and the sub-peak electron content, re- 
spectively. These parameters have been displayed 
for separate months in the latitude-local time plane, 
in the NBS Technical Note Series 40; examples of 
them have been discussed above. 

With one year of data available, it becomes prac- 
tical to display this large quantity of information in 
very compact form. For each of these four param- 
eters, and each of the five stations, we have prepared 
contour ' 'maps' ^ showing the diurnal versus seasonal 
variation of these parameters. Figures 8 through 11 
may be referred to in the discussion which follows. 

We shall first discuss the various parameters 
individually from a phenomenological point of view, 
drawing attention to their characteristic diurnal 
and seasonal features and to the latitude variation 
evident from the five stations' data. A synthesis 
and interpretation is then attempted in section 5. 

4.1. Height of the F2 Peak {hmaxF2, fig. 8) 

Characteristically, hmeixF2 is lowest at the time 
immediately following sunrise, as low as 250 km in 
the winter and at higher latitudes. It rises through- 
out the daylight period and generally continues to 
rise after sunset and throughout the nighttime period. 
The total diurnal excursion of Amaxi^2 is surpris- 
ingly constant (100 km), whatever the season or 
latitude, except that at midlatitudes the variation is 
somewhat less in winter. There is a distinct similar- 
ity between the daytime 300 km contours at the 
station pairs White Sands/Grand Bahama, and Fort 
Monmouth/Newfoundland. These pairs of stations 
have nearly the same magnetic dip (60° and 72°, 
respectively) , although their geographic latitudes are 
rather dissimilar. 

The largest values of hnvcixF2 occur in the summer 
nighttime, and are typically 400 km. At all but the 
lowest latitude studied, these high values are centered 
on local midnight; at Puerto Rico, they occur before 
midnight, centered on 2100 local time in the summer 
months. It is interesting to note that this tendency 
is entirely absent at Grand Bahama, which is farther 
north by only 8° (geographic, geomagnetic, or dip). 
On the other hand, the daytime similarity noted 
above between stations with similar magnetic dip, 
is also evident in the diurnal/seasonal patterns of 
the nighttime 400 km contour. 

4.2. Quarter-Thickness of the F2 Peak {Scat, fig. 9) 

The possibility of interpreting this quantity in 
terms of the (neutral particle) scale height at the F2 
peak (see section 2.3) suggests that Scat may be an 
important indicator of the temperature in the F 
region. We shall keep this possibility in mind, as 
the phenomenology of Scat is discussed. 
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Figure 8. Diurnal and seasonal variations of hmax. 
Contours in km. 



Scat shows a diurnal variation at all latiludos and 
seasons, and this is most marked in the high latitude 
summer. Noon values of Scat approach 100 km in 
the high latitude summer, but the niglittime values 
are more typically 50 km or less. At Newfoundland, 
the hig:hest latitude, the sense of tlie diurnal variation 
is reversed in the winter: maximum values then occur 
in the pre-dawn hours, and generally occur earlier 
in the equinox months than at the winter solstice. 
At all ot her latitudes, the daytime values are distinctly 



larger than at night. As w^ith AmaxF2 there is 
somewhat more similarity between stations of similar 
magnetic dip, than among those of similar geographic 
latitude. 

There is a smaller seasonal variation at the lowest 
latitude (Puerto Rico), with minimum values in the 
winter and maxinnmi in the sunnner. The total 
change is less than at high latitudes, but the smallest 
values (30 km, winter pre-dawn) are as small as the 
minima elsewhere. 
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SCAT NEWFOUNDLAND 




FORT MONMOUTH 
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GRAND BAHAMA 



PUERTO RICO 





Figure 9. Diurnal and seasonal variations of Scat. 
Contours in km. 



There is a suggestion of a semi-diurnal variation in 
midlatitude winter, with maxima at 0200 and 1400, 
minima at 0730 and 1930. 

4.3. F2 Maximum Electron Density (iVmaxF2, fig. 10) 

Although the phenomenology of this parameter is 



well-laiown from its direct relationship with foF2, 
we shall review its phenomenology here in order to 
complete our picture of variations at the F2 peak. 
It will also prove interesting to compare A'maxi^2 
with the F2 subpeak electron content. The very 
great range of iVmaxF2 (through a factor of ten at 
all latitudes) renders the maxima very prominent, 
and the minima only slightly less so. 
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N max NEWFOUNDLAND 
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GRAND BAHAMA 



PUERTO RICO 





Figure 10. Diurnal and seasonal variations of NmaxF2. 
Contours in (electrons/cm 3) XI 0"^. 



In the daytime, high latitudes show a single winter 
post-noon maximum, which shifts wildly to late 
afternoon in summer, while at lower latitudes the 
post-noon maximum is between 1300-1400, nearly 
independent of latitude or season. The similarity 
between the late afternoon summer maxima at Ft. 
Monmouth and Newfoundland suggests magnetic 
rather than solar control for this anomaly. 



Low latitudes show distinct equinoctial maxima, 
with minima in the smnmer and winter. At inter- 
mediate latitudes, the daytime seasonal variation is 
complicated by the transition between these two 
clear patterns. The equinox noon maxima at low 
latitude provide the very highest values of iVmax. 
At higher latitudes, the highest values are nearly 
independent of latitude, but possess the diurnal and 
seasonal pattern just mentioned. 
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Figure 11. Diurnal and seasonal variations of Shmax. 
Contours in (electrons/cm 2 col)X10-i2. 



At night, the minimum values always occur in 
winter, and the largest values occur in summer, at 
all latitudes. The minimum value of A^max is 
nearly the same at all latitudes, and occurs in the 
pre-dawn hours. The period of time occupied by 
the minimum (seasonally centered on December and 
diurnally centered on about 0500), becomes larger 
with increasing latitude. 



4.4. Subpeak Electron Content (S/imax, fig. 11) 



The detailed variations of Shm^x show the in- 
fluence of each of the preceding parameters. Again 
considering first the daytime periods, it is clear that 
the pattern of the maxima is the same as for iVmax: 
high latitudes display a winter maximum which 
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''divides'^ into two equinoctial inaxima at lower 
latitudes. Tlic maxima are gcMUM'ally post noon, but 
not as late as those of A^max. The pre-noon increase 
in Skmiix is more uniform with season and with 
latitude than is the case with iVinax. The maxi- 
mum values of Shnuix decrease witii increasing 
latitudes. 

At night, the minimum values of SAmax occur in 
the pre-dawn period, and are lowest in winter. 
However, there is a tendency, also suggested in the 
A^max data, for the lowest values to occur before 
the winter solstice, rather than at it. The minimum 
values of Shmo^x decrease with increasing latitude. 



5. Normal Conditions in the F Region; the 
Seasonal Anomaly 



An understanding or interpretation of the com- 
plicated behavior of the /"region, as illustrated above, 
is handicapped from the start by the fact that no 
shigle influence appc^ars to exert dominance in a 
consistent way. While solar radiation is obviously 
the ultimate influence, it is well known that the 
elementary Chapman theory of ionospheric layer 
formation predicts that the maximum electron den- 
sity should vary as (cos x)^, and it is equally well 
known that the F2 layer does not obey this law. 
However, Eatclifl'e [1951] observed that somewhat 



better correspondence with this law was obtained if 
one used, instead of A^max, the subpeak electron 
content ShmsiX. 

Seasonally, one of the anomalies of the F2 layer 
is that iVmaxi^2 is greater in winter than in summer, 
quite in contradiction of simple theory. This has 
been illustrated already in the discussion of the 
diurnal and seasonal patterns of iVmaxi^2 (section 
4.3). "Ratcliffe^s [1951] data suggested that this 
anomaly was not present in SAmax, and his con- 
clusion was supported (but imperfectly) by Osborne 
[1952]. We should now like to demonstrate that 
the cos X variation of Shimix is indeed more regular 
than that of A^max, but that a seasonal anomaly 
persists in Shmux hi mid and higher latitudes. 
Furthermore, our data show that it is the summer 
which should be considered anomalous. 

Consider the lower portions of figure 12. The 
mean quiet values of A/max between the hours of 
sunrise and noon are plotted logarithmically versus 
cos X for the year's data, May 1959-April i960, for 
the two stations Newfoundland and Puerto Ilico. 
Summer values are shown by X's, winter values by 
open circles, and equinox values by filled circles. 
The seasonal anomaly mentioned above is evident 
in the considerable extent to which the winter values 
of A/max exceed those in summer. The anomaly is 
less marked at Puerto Kico than at Newfoundland. 

On the other hand, consider the upper part of this 
figure. Here, the corresponding variations of SAmax 
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Figure 12. Variation of Shmax {upper portion) and Nmax (lower portion) versus cos X) 

between hours of sunrise and noon at Newfoundland and Puerto Rico. 

Key: winter values Oi equinox values •, summer values x; no distinction made for Puerto Rico ShmsiX. 
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versus cos x are shown for the same stations. It is 
at once apparent that the seasonal anomaly is 
entirely absent at Puerto Rico, but easily detectable 
at Newfoundland. At the higher latitude it is, 
however, distinctly less well marked than in A'^max. 

Stations at intermediate latitudes (Grand Bahama, 
White Sands, and Fort Monmouth) were examined 
in a similar way, and found to behave in a manner 
intermediate between the extremes shown here, in 
correspondence with latitude. 

The line passing through the Puerto Rico ShmsiX 
versus cos X data is also drawn on the corresponding 
graph for Newfoundland, and it is seen that it agrees 
reasonably well with the winter points. This same 
agreement is found with the winter data at the inter- 
mediate stations mentioned above. Therefore, it 
appears that it is the summer season which is anom- 
alous, and that in the winter time the subpeak 
content of the ionosphere approximately obeys the 
Chapman Law at all latitudes. It is remarkable that 
at the latitude of Puerto Rico, this same law is 
obeyed at all seasons. The line in this figure is 
given by SAmax^ 2.2X10^^ (cos x)^-^^ electrons per 
cm^ column. In this one regularity of the subpeak 
electron content, it appears that we have an impor- 
tant clue to some of the causes of the otherwise 
puzzling behavior of the F region. 

6. A Synthesis and Partial Interpretation of 
the Preceding Observations 

6.1. Assumptions 

It would be quite difficult, and perhaps impossible 
at this point, to give a quantitative explanation for 
all the features of the mean quiet ionosphere por- 
trayed in the preceding figures. However, a quali- 
tative discussion is certainly useful, since the various 
factors to be mentioned must at least account quali- 
tatively^ for these data before quantitative tests can 
be devised. We shall therefore attempt to construct 
an argument which can account qualitatively for 
most of the behavior described above. At the outset, 
several assumed conditions and relationships should 
be stated: 

(1) It is assumed that the main features of figures 
8-11 are representative of the real ionosphere, and 
are not seriously influenced by errors in the calcula- 
tion of the individual N{h) profiles. 

(2) The calculated F2 layer quarter-thicknesses 
(Scat) are taken to indicate the variations in neutral 
particle scale height at the level of the F2 peak as 
discussed in section 2.3. 

(3) Thermal equilibrium is assumed, and therefore 
the variations in Scat are taken as indications of 
the temperature variations of the neutral atomic and 
molecular species at the level of the F2 peak. 

(4) The height of the F2 peak is presumed to be 
influenced most strongly by diffusion, balanced by 
an attachment-like loss process [Ratcliffe, 1960, 
p. 440]. 

(5) The attachment-like loss process has an 
effective attachment coefficient assumed to be of 



the form 13=^0 exp { — (h—ho)/Hxy], where H^y is the 
scale height of the molecular species involved in the 
ion-molecular charge exchange process, probably 
molecular oxygen [Ratcliffe, 1960, p. 383]. 

6.2. Diurnal and Seasonal Temperature Variations 
at the F2 Peak 

By assumption 3 above and the diurnal/seasonal 
variations of F2 quarter-thickness discussed in 
section 4.2, it is evident that large diurnal, seasonal, 
and latitudinal temperature variations occur at the 
level of the F2 peak. It is possible that a small 
part of this variation may be ascribed to vertical 
motion of the layer in a region with a height gradient 
of scale height, but it seems clear that large tem- 
perature variations must occur throughout the F2 
layer if the layer thickness is an}- guide. Other 
evidence for the diurnal variation has recently been 
given by Jacchia [1961]; he has described diurnal 
variations of F region temperature as determined 
from satellite drag data that are similar to the 
variations suggested here. 

However, the scale height, or temperature varia- 
tions implied by figure 9 are not attributable solely 
to solar radiation as the heat source. In fact, the 
following evidence of two significant sources of heat 
at the F2 peak may be seen in figure 9. 

(a) Seasonally, diurnally, and latitudinally, the 
temperature increases with decreasing solar zenith 
angle, consistent with solar heating of the atmos- 
phere. 

(b) Data from the higher latitude stations suggests 
that corpuscular heating is becoming important 
above about 70° dip. Evidence for this is the 
following : 

i. Above 70° dip, the summer daytime F2 
thickness is at least 25 percent larger than in 
similar periods at lower latitudes, despite the 
smaller zenith angles there. This difference can- 
not be explained by a gradient of scale height with 
height, since hmeixF2 is slightly lower at higher 
latitudes. 

ii. The minimum nighttime F2 thickness at 
the highest latitude is also about 25 percent 
greater than the corresponding values at lower 
latitudes. If this were to be explained entirely by 
the slight difference in nighttime Amaxi^2 between 
the two latitudes, it would require an improb- 
ably large height gradient of scale height of 
dHldh= 0.5. 

iii. In high latitude winter, the F2 thickness is 
larger at night than in the daytime, reaching its 
diurnal peak in the early morning hours, and at 
an earlier hour in the equinox. This behavior 
corresponds to the diurnal and seasonal incidence 
of auroral activity at these latitudes [Harang, 
1951]. 

Thus, these data suggest the following general 
thermal behavior: at low latitudes the heating of the 
F2 region is predominantly solar, and depends pri- 
marily on the solar zenith angle. At higher latitudes 
there is superimposed on the solar heating, an auroral- 
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zone heating of the F2 region that increases all 
teniperatures by roughly 25 percent. 

6.3. Effect on the Production and Loss of Ionization 

The increases in temperature described above 
cause increases in the neutral particle density at F2 
heights. In an atmosphere in diffusive equilibrium 
containing atoms and molecules of the species and 
O2, an increase in the total density implies an increase 
in tlie atomic density that goes as the square root 
of the increase in the molecular density, if the rate 
of production of the former from the latter is not 
altered.^ 

Since it is believed that the rate of loss in the F 
region varies as 71(02) while the rate of production 
varies as 7i(0), it can be seen that an increase of 
temperature increases the loss rate to a greater 
extent than it permits an increase in the rate of 
electron production. 

6.4. A Qualitative Explanation of the Seasonal 

Anomaly 

a. Daytime 

At low latitudes it appears that the electron pro- 
duction and loss processes, even though influenced 
by dim^nally varying temperature, cause a '^ normal' ' 
electron density variation that is approximately 
Chapman-like in the morning hours at all seasons, 
although the reasons for this remain to be explained 
quantitatively. At higher latitudes in the summer- 
time, auroral and solar heating together increase the 
loss rate over the winter value to such an extent tliat 
the overall electron content in daytime is signifi- 
cantly reduced, despite the fact that the electron 
production rate has increased also. During high 
latitude winter, the solar heating component is 
reduced and the remaining auroral heating compo- 
nent is apparently insufficient to destroy completely 
the ^'normar^ behavior of the F region. It should 
be noted that it is the summer season during which 
the Fl layer appears. 

b. Nighttime 

At night, there is no seasonal anomaly: Night- 
time densities are smaller in the winter than in the 
summer and this is true at all latitudes, as may be 
seen in the iVmax and SAmax diagrams of figures 
10 and 11. At first sight, this may seem inconsistent 
with the explanation just advanced for the daytime 
variations. However, consideration of the factors 
dominant at night shows that this is not so. The 
nighttime-variations of electron content are explain- 
able in terms of the equilibrium height of the F 
region affected by combined processes of diffusion 
and attachment-like electron loss [Duncan, 1956]. 
By this theory, the nighttime F region peak settles 
rapidl}' to an equilibrium height at which the effec- 
tive attachment coefficient is given by 



&=g 



sm"^ (p 
2^ 



where g is the gravitational acceleration, H is the 
scale height of the atomic species at the level in 
question, v is the positive ion collisional frequency, 
and ip is the magnetic dip. We shall call the value 
of /3 at this height the ^'equilibrium value,'^ jSe. 
At the same time, the electron density distribution 
approaches the Chapman form 



where 



A'^A^max exp ^{l—z—e~^) 

2= r dh/H. 

J Jimax 



(4) 



(3) 



2 1 am indebted to my colleague R. B. Norton for emphasizing the importance 
of this. 



These same arguments justify the interpretation of 
Scat in terms of the scale height H, and provide the 
basis upon which the extrapolations above Amax F2 
have been made. 

Equation (3) indicates that an increase in scale 
height, or temperature, causes a decrease in the 
equilibrium value of /3e. The positive ion collisional 
frequency, which varies directl}^ with the neutral 
particle density, will also increase contributing to the 
diminution of the equihhrium value, jSe.) At the 
same time, the values of 13 existing in the atmosphere 
are larger iFor larger H, as explained in paragraph 6.3 
above. Two effects therefore follow from the higher 
summer temperatures: (a) the nighttime F region 
will be found at a greater equilibrium height such 
that /? has its required smaller value; and (b) the 
summer nighttime F region will decay at a slow 
rate corresponding to this smaller value of 13, result- 
ing in tlie maintenance of higher electron densities 
in summer than in winter. In the Amax diagrams 
of figure 8, the nighttime summer equilibrium 
heights (at, say, midnight) are seen to be higher than 
the winter lieights, in agreement with this argument. 

Referring again to the diurnal/seasonal variations 
of Shmux in figure 11, it can be seen that the values 
of Shnvdx at sunset do not vary greatly with season 
at any latitude. The faster winter loss rates at the 
lower equilibrium heights will then rapidly decrease 
the electron content (and Mnax), thereby accounting 
for the lower values of these parameters in winter. 

7. Summary 

Our examination of one 3^ear of mean quiet electron 
density data between latitudes of 18° N and 47° N 
has been concerned primarily with the properties of 
the i^ region near its peak, together with the subpeak 
electron content of the layer. We have found that, 
if the characteristic thickness of the region near the 
peak is interpreted in terms of the local scale height, 
and hence in terms of the local temperature, the 
variations of height, thickness and electron content 
may be explained qualitatively. Several well known 
anomalies of the F region are similarly accounted 
for — the seasonal, diurnal, and latitudinal anomalies. 
The ^Vinter'' anomaly, in particular, is found to be 
misnamed; the data suggest strongly that it is the 
summer period which should be considered anomalous 
for its low electron content. The seasonal temper- 
ature variations together with an auroral-zone 
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heating component suggested by the data, are used 
to explain this summer anomaly in terms of an in- 
creased electron loss rate; it is suggested that the 
Fl layer is similarly explainable. The absence of a 
summer anomaly at night is found to be explainable 
by the loss rate at the height of the layer peak: 
with an increase in temperature the layer seeks a 
greater height where the loss rate is smaller. 



As implied in the introduction, this study has been 
possible only because of the assistance rendered at 
the data's primary source, the ionospheric field 
stations; I wish to express my sincere appreciation 
for the efforts of the many station scientists involved. 
The NBS N(h) group led by G. H. Stonehocker, has 
rendered a service that is second only in chronology. 
Similarly invaluable has been the programing of the 
work for NBS' high speed computer by Dr. H. H. 
Howe. Discussions with the author's colleagues 
S. Kadicella, T. N. Gautier, R. A. Duncan, and R. B. 
Norton have been most helpful in delineating the 
ideas offered here. 

8. Appendix. A Note on the Day-to-Day 
Variability of the Ionosphere 

As was mentioned earlier in the discussion of the 
averaging process, the standard deviations (a) and 
relative standard deviations (o-/mean) of the profile 
data are also obtained. These parameters give a 
quantitative measure of the day-to-day variability 
of the quiet ionosphere, and permit some interesting 
conclusions regarding the degree to which various 
external influences control the ionosphere. 

The standard deviation of the (roughly 20) 
values of Amax, Shmax, Scat, and A^max entering 
their respective means have been determined hourly 
for each month of the year May 1959-April 1960. 
In the case of the electron density parameters 
Shmhx and iVhiax, it is immediately evident that 
the standard deviations are closely proportional to 
the mean values themselves, and that a clearer 
picture of the variability of the quantity is obtained 
if the standard deviation is expressed in percent of 
the mean value. This is also meaningful for the 
quarter-thickness. Scat, but has little meaning for 
Amax. Accordingly, our figures and remarks will 
concern the percentage variability of ShuiSiX, A^max 
and Scat, but will deal with the true variability of 
hmsix, in kilometers. 

From the one-year's data available to this study, 
it has not proven possible to observe any important 
seasonal variation in these variability data. There 
is a slight tendency for winter daytime values to be 
smaller (especially for A^max and Hma,x) than at 
other times and seasons, but this has not been 
explored further. It is interesting to note, however, 
the conclusions of section 5 in which other evidence 
suggests that winter shows a generally less anomalous 
behavior than summer. 



Because of the relatively minor seasonal depend- 
ence of the variability parameters, all months of the 
year's data have been averaged together for the 
preparation of the latitude versus local time maps 
of figures 13, 14, 16, 17. These maps show contours 
of the true variability of hmax (in km) and the 
percent variability of Scat, Shmsix, and iVmax, 
respectively, between latitudes of 15° N to 50° N, 
versus 75° W time. 

8. 1 . Variability of hmax 

As shown in figure 13, Amax F2 has about twice 
the variability at night as in the daytime, and is 
generally larger at lower latitudes. The minimum 
variability is at midday, but the maximum variability 
is generally post-midnight — markedly so, at lower 
latitudes. 

8.2. Relative Variability of Scat 

The percent variability of the F2 layer quarter- 
thickness is shown in figure 14. Scat has been inter- 
preted as a measure of the neutral scale height and 
hence the temperature at the F2 peak throughout this 
paper. Its variability is generally larger at high 
latitudes and in the early morning or sunrise periods, 
ranging from 10 percent to 30 percent of its mean 
value. The greater variability at high latitudes may 
be related to the importance of corpuscular heating 
at latitudes approaching the auroral zones, as pro- 
posed in section 6.2; it is reasonable that any corpus- 
cular heating of the atmosphere will have more day- 
to-day variability than the heating due to solar 
radiation. 

It is of interest to inquire to what extent the 
variabihties of the height of the F2 peak and its 
characteristic thickness are related. It might be 
imagined, for example, that all the variability of 
hmhx might be accounted for by variability in the 
thickness of the layer. This would be the case if the 
layer were essentially fixed in position and merely 
varied in thickness from day to day. Alternatively, 
it might be imagined that the whole region ^^moves^' 
vertically from day to day, producing a variation in 
hmsix that is essentially independent of variations in 
the layer thickness. 

To examine this point, in figure 15 we have plotted 
diurnally the standard deviations (in km) of /imax 
F2 and Scat for the entire year's data and for the two 
stations Newfoundland and Puerto Kico. The 
shaded area shows in each case the part of the vari- 
ability of Amax which cannot be accounted for di- 
rectly by the variability in Scat, or vice versa. 

First, one notes that generally Amax varies more 
than Scat, except during high latitude daytime, when 
the two are essentially the same. There is a great 
difference between the two at night, indicating that 
the layer is found at considerably different pX)sitions 
(altitudes) from day to day with comparatively 
little change in thickness. There is, however, some 
indication of a correlation between the magnitudes 
of the two parameters. 
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DIURNAL, LATITUDINAL VARIABILITY OF HEIGHT 
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Figure 13. Diurnal, latitudinal variability of hmaocF2; quiet 
7nean values averaged over May 195 9- April 1960; contours 
in kilometers. 
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Figure 14. Diurnal, latitudinal relative variability of Scat; 
quiet mean values averaged over May 195 9- April 1960; 
contours in percent. 
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Figure 15. Comparative variability of hmaxF2 and Scat 
versus time, at Newfoundland and Puerto Rico, averaged over 
May 1959-April 1960. 



8.3. Relative Variability of Total Electron Content 

This parameter (fig. 16) has a very regular 
seasonal and time dependence. It is nearly constant 
(30%) at night, independent of latitude, and 
appears to depend ver}^ simply upon the solar 
zenith angle, being smallest chu'ing daytime at low 
latitude. The diurnal minimum vahu* occurs, 
however, from one to three hours postiioon. 

8.4. Relative Variability of Maximum Electron 
Density 

In contrast with the variability of Shnmx, the 
variability of the maximum electj'on density (fig. 17) 
suggests a mixture of latitudinal and solar control. 
The maximum values (40%) occur in the predawn 
period at high latitudes. The minimum values 
occur towards lower latitudes in the daytime, and 
these become equal to the relative variability in 
subpeak electron content discussed in paragraph 3. 
This is consistent with the Chapman-like behavior 
of the low latitude daytime ionosphere discussed in 
section 5 : If the F region possesses a Chapman-like 
behavior, the electron density at all altitudes varies 
in proportion to the maximum density, and hence 
the electron content does also. This behavior is 
somewhat evident at all latitudes near midday, but is 
definitely not the case at night — especially at higher 
latitudes. However, we may note that the high 
latitude predawn peak variability of ShiasiX occurs 
at nearly the same time as the peak variability in 
layer thickness (vide supra). The behavior of the 
variabilities of Mnax, SAmax, and Scat at high 
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DIURNAL, LATITUDINAL RELATIVE VARIABILITY OF 
SUB PEAK ELECTRON CONTENT 
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Figure 16. Diurnal, latitudinal relative variability of Shmax; 
quiet mean values averaged over May 195 9- April 1960; 
contours in percent. 



DIURNAL, LATITUDINAL RELATIVE VARIABILITY OF 
F2 MAXIMUM ELECTRON DENSITY 



QUIET DAYS, MAY 1959-APRIL I960 




Figure 17, Diurnaly latitudinal relative va? lability of NmaxF2; 
quiet mean values averaged over May 195 9- April 1960; 
contours in percent. 



latitudes during nighttime are again consistent witli 
the corpuscular heating discussed in section 5.2: 
As the day-to-day changes in temperature modulate 
the thickness of the nighttime jF region, its maximum 
density varies in inverse fashion without appreciable 
change in the electron content. 

8.5. Summary 

The day-to-day variability of an ionospheric 
parameter, expressed as a percentage of the monthly 
mean, has been shown to be subject to interpretation 
consistent with the physical processes controlling 
the layer. The variabihty of a given parameter 
shows properties relatable to certain physical 
influences which may be obscured in the real values 
of the parameters themselves. In particular, the 
variability of the high latitude nighttime ionosphere 
is consistent with a variable corpuscular influence 
(probably heating) not evident at lower latitudes. 
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